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, Abstract. We analyze the physical justification of the picture proposed 

fSj ' by Salvador-Sole et al. in these proceedings for the time evolution of the 

universal density profile of dark-matter halos. According to this picture, 
halos have at formation a stable (i.e. independent of mass and time) di- 
mensionless density profile, the characteristic length and density scales 
| of the profile depending on the underlying cosmogony. Subsequent evo- 

■"sj- . lution is driven by mass accretion onto the outskirts of halos and can be 

characterized simply by the increment of halo radius with time and the 
corresponding decrease of the critical density of the universe. We find 
| this picture to be a reasonable good description of the expected evolution 

■ of halos in hierarchical models of structure formation. 
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1. Introduction 



Relying on data from high resolution N-body simulations performed by Navarro, 
Frenk & White (1997, NFW), Salvador-Sole et al. (1998, S98) derived the laws 
^ ■ for the evolution of the dimensionless characteristic length and density scales 

(x s and 5 C respectively) of the universal density profile of dark-matter halos 
proposed by NFW. For a halo of mass M at time t (redshift z), those laws have 
the following expressions: 

x s (t,M) = x sf ^, (1) 

5c ^ M) = 6ci n( Zi )(i + z)> (2) 

where R(t) is the virial halo radius at t, O(^) is the cosmic density at z in units 
of the critical density of the universe (p C rit)j and if and Zf are the formation time 
and formation redshift of the halo, respectively. It is important to remark that 
equations ([l]) and (^) are two independent fitting formulae. N-body simulations 
show that x s and <5 C are, in fact, linked through the condition that the mean 
internal density of halos within the virial radius is equal to about 200p cr it (im- 
plying that halo density profiles are one-parametric functions). Equations (ffl) 
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and (g) are found to be consistent with that property for different cosmological 
models (cf. S98). 

It is evident that the proportionality constants x s f and <5 c f (which are cos- 
mogony dependent) correspond to the values of x s and 8 C when halos form, 
respectively. Therefore, the values of the scale radius, r s = x s R(t), and of the 
characteristic density, p c = S c p C n t (z), of halo density profiles are set at the halo 
formation time. Halo formation is basically characterized by the last major 
merger yielding a substantial re-arrangement of the internal structure of merg- 
ing halos (see Salvador-Sole et al. in these proceedings for a quantitative defi- 
nition). After formation, subsequent evolution by matter accretion (secondary 
infall and/or minor mergers) does not change the values of r s and p c . However, 
since pcrit(^) decreases with time while R(t) increases accordingly, the values of 
x s and 5 C change as described by equations and (Q). 

In the present work we develop a physical model with the aim of checking 
the validity of the picture proposed for the evolution of the scaling parameters 
of halo density profiles. On the other hand, we will try to determine for a set 
of cosmogonies analyzed the corresponding shape of the stable density profile of 
halos at their formation time. 

2. The Structure of Halos Formed by Binary Major Mergers 

In this section we describe the model linking the scaling parameters of the density 
profile of new formed halos and their progenitors. Let us assume that halos form 
essentially through the major merger of two halo progenitors (this approach is 
justified in the next section). In practice, as a major merger we mean that 
if a halo has a formation mass Mq, then its two halo progenitors have masses 
Mi and M 2 (with Mi > M 2 ) so that M 2 /Mi > A m = 0.6. The threshold A m 
between major and minor mergers is an empirical parameter introduced to allow 
for a better motivated definition of the formation time of halos. The value of 
A ffl = 0.6 has been derived by fitting the mass-density correlation of halos in 
high resolution N-body simulations (see S98). 

Consider a halo of mass Mq with internal energy Uq at its formation time, 
if, and the system formed by its two halo progenitors at the time they reach 
turnaround. The total energy of the system at turnaround can be written as: 

E ta = U X + U 2 + E 12 , (3) 

where f/j is equal to the sum of the internal energy of the virialized mass of 
the i-th halo progenitor at turnaround plus the total energy of the mass that 
will be accreted onto that halo from turnaround until if. The last term in the 
right-hand-side of equation @ gives the mutual gravitational interaction plus 
orbital kinetic energy of the system at turnaround. 

If there is not significant mass-loss neither during the formation process of 
Md, nor during mass accretion onto halo progenitors from turnaround until if, 
then from energy conservation one has 

Uq = Ui + U 2 + E 12 , (4) 
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and also for i = 1,2 one can approximate f7j by the internal energy of halo 
progenitors at if . Therefore, for i = 0, 1, 2 one can write 



where Mj and Ri are the mass and the radius of each halo at if, respectively, 
and Mi + M2 = Mq. In writing equation (|5|) it has been assumed that halos are 
virialized systems described by a universal density profile, which determines the 
specific expression of function F(x s ). 

Concerning the expression of E12, it can be written as 
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where D m is the turnaround separation between halo-progenitor centers. In the 
well known point-mass approximation and taking if equal to one orbital period 
one gets: 

Dl = (2-SfGM ^y . (7) 

The so-called circularity parameter, S, appearing in equations (^) and (0) takes 
into account that halo progenitors can merge following a non-radial orbit. For 
an elliptical orbit of eccentricity e, one has S = 1 — e. Thus, a radial orbit 
corresponds to S = while a circular orbit to S = 1. For an object moving in 
an arbitrary orbit of total energy E one has S = J 2 /Jc( E ) (Merritt 1985), where 
J is the angular momentum of the orbit and J C (E) is the angular momentum 
of a circular orbit of energy E. Non-radial mergers are caused by external 
torques from the large-scale matter distribution surrounding halo progenitors. 
In hierarchical structure formation scenarios one expects that low-mass halos 
will merge in non-radial orbits while massive halos will tend to merge in nearly 
radial orbits. On the other hand, the mutual tidal interaction between merging 
halos slows-down the orbital motion and the final merger time will be longer 
than the value computed for point masses (eq. 0]). In practice, this can be 
accounted for by assuming a non-null value of S, which for a fixed D m increases 
if by a certain amount. The larger the mass of merging halos the stronger the 
expected effect of mutual tides on the merger time. Therefore, the combined 
action of external torques and mutual tides between halo progenitors can be 
described by a value of S different from zero for a wide mass range. The exact 
mass dependence of S is difficult to predict, but in a first approximation the 
influence of both external torques and mutual tides on the final results can be 
addressed by taking a constant value of S over the whole mass range analyzed. 

Substituting equations (||), (]6|), and (^) into equation (Q) one derives the 
following equation: 



where r = t{H(t{) and H(t) is the Hubble parameter at t. In the left-hand-side of 
equation (§) we have written x s { instead of x s q because the latter corresponds to 
the value of the dimensionless scale radius of Mq at its formation time. The ex- 
pression of F(x s ) for spherically symmetric halos in hydrostatic equilibrium with 
an isotropic velocity distribution is completely determined by the halo density 
profile, p(r), which is assumed universal. The different shapes for that univer- 
sal profile considered in the present work are particular cases of the following 
general expression (Zhao 1996) 



p(r) 



ya(l + y/3)V ' 
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where y = r/r s and r is the radial distance to the halo center. From this density 
profile one can infer the corresponding mass profile M(r), the one-dimensional 
velocity dispersion profile cr(r), and finally, one can compute F(x s ) from the 
following expression: 



F(x s 



U 



GM 2 /(2R) [m(c)p 



m{y)p(y)ydx - c 3 p(c)a 2 (c) 
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where p(y) = p(r)/p c , m(y) = M(r)/(47rr s 3 /j c ), &(y) = a(r)/y/4irGr^p c , and 
c = x" 1 is the halo concentration parameter. 

To sum up, we detail the steps followed by the practical implementation of 
the present model: 

Step 1: Starting from a halo of mass M at a fixed time t, we compute by 
means of the S98 clustering model the halo formation time tf , its formation mass 
Mq, the typical mass of its halo progenitors, M\ and M2, and the formation time 
of the latter, tfi and tf2- 

Step 2: An empirical value of the dimensionless scale radius (denoted by 
Xg f ) is assumed for halo progenitors of any mass at their formation time. This 
is equivalent to setting a stable density profile at formation for all progenitors 
considered if a universal halo density profile is assumed. 

Step 3: The value of is evolved from the formation time of each halo 
progenitor ta till their merger time corresponding to the formation time tf of 
the resulting new halo Mq. This yields the dimensionless scale radii for the 
density profile of halo progenitors, x s i and x S 2, at time if that will be inserted 
into equation (|). The evolution of x^ f can be computed from either equation 
(|l]) or equation (2|). We have checked that both yield consistent results to within 
typically 10%. 

Step 4: Solving equation (@) for x s f we derive the theoretical dimensionless 
scale radius for the density profile of Mq at formation (denoted by x* f ) which 
has to be compared with the corresponding empirical value assigned to halo 
progenitors in step 2, x e si . 



3. Progenitor Masses and Binary Major Mergers 

As stated at the beginning of the previous section, we assume that major mergers 
tracing the formation of new halos are essentially binary. This can be justified 
in the framework of the clustering model developed by S98. 
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Consider a halo of mass Mi that incorporates a certain amount of mass 
AM at time t. According to S98, if AM /Mi > A m then Mi is destroyed and 
a new halo of mass Mq = Mi + AM is formed by a major merger, while if 
AM /Mi < A m the halo M = M x + AM is the result of the evolution of the 
halo Mi by mass accretion and/or minor mergers, hence Mi is not destroyed, it 
simply changes its mass to Mq. 

Now assume that, in the case of major mergers, the mass AM always comes 
from a single halo that hereafter we denote by M%. It is easy to see that 

A m M Mq 

t — < Mi < — , for i = l,2. (11) 

1 + A m 1 + A m ' K J 

This interval limiting the mass of halo progenitors is centered at Mq/2, and 
| Mi — Mq/2\ = \M2 — Mq/2 I , or equivalently, Mi and M2 are at the same 
distance of Mq/2. Therefore, the distribution function of halo progenitor masses 
should be a symmetric function around Mq/2. This distribution function can be 
derived from the S98 clustering model and, hence, its symmetry can be checked 
(see Raig et al. 1998). We find that indeed the symmetry is fulfilled to within 
3% and this result is independent of redshift and of the cosmogony considered. 
The symmetry of the distribution function of halo progenitors strongly supports 
that major mergers are essentially binary. 



4. Results 

In this section we will present the results of the comparison between the dimen- 
sionless scale radius at formation for halos of mass M and their progenitors (x* f 
and respectively), according to the model described in §2. This comparison 
will be performed for different cosmogonies and a fixed universal density law 
of the form given by NFW. Subsequently, we will consider the implications of 
assuming a different universal density profile. The set of cosmogonies analyzed 
in the present work is described in columns (1) to (5) of Table 1. We have 
considered a standard biased cold dark matter model (SCDM), a flat cold dark 
matter model with non-null cosmological constant (ACDM), and a power-law 
model with power index n = — 1 and flat geometry. These models were also 
analyzed by NFW and S98. In Table 1, Qq is the matter density parameter, Ao 
is the cosmological constant in units of 3Hq, where Hq = 100 h km s _1 Mpc -1 , 
and <7g is the rms mass fluctuation in spheres of radius 8 Mpc. 

In Figures |l] and ^ we show the fractional difference between x* f and 
for the SCDM and the ACDM models, respectively, assuming a universal NFW 
density profile. Similar curves are found for the power-law model. Results 
are plotted as a function of halo masses at z = scaled to the characteristic 
mass M* defining a unity rms density fluctuation. For cosmogonies which are 
self-similar or close to self-similar the results expressed in this way are redshift 
independent. Thick solid lines correspond to different values of the parameter 
S and were derived assuming a value of _c|- equal to the best fit dimensionless 
scale radius found by S98 from their analysis of the NFW N-body simulations. 
The concrete values are shown in column (7) of Table 1, and in column (8) are 
listed the 90% confidence intervals resulting from the x-square fit. 
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Figure 1. Fractional difference between the dimensionless scale ra- 
dius at formation for halos of mass M and their progenitors (x l si and 
Xg f respectively) in the case of an SCDM cosmology and assuming an 
NFW universal density law for dark-matter halos. The different thick 
solid lines correspond to different values of the parameter S ranging 
from S = (upper curve) to S = 1 (lower curve) in steps of 0.25. They 
were derived by taking equal to the best fit value obtained in S98 
(see Table 1). Thin solid lines correspond to twice (upper) and half 
(lower) the latter value for S = 0. Finally, the thick dashed line has 
been derived for the value of Xg f that minimizes the maximum frac- 
tional difference in the whole mass range (denoted as x 1 ^ in the text). 
The shape of this latter curve is independent of the value of S assumed. 
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Figure 2. Same as Fig. |T| but for the ACDM model. 



As can be seen from Figures [l] and |2| the fractional difference is not far 
from zero as expected if the dimensionless scale radius of halos were nearly fixed 
at formation. It is important to notice that taking into account the dynami- 
cal effects of external torques and mutual tides (S > 0) tends to improve the 
consistency with a universal value of x B {. The maximum fractional difference is 
typically between 20% and 35% and is found for low-mass halos. For massive 
halos the agreement is quite satisfactory. On the other hand, another interesting 
result is that our model can be used to predict the value of which is in better 
agreement with the stability of halo density profiles at formation. As we show 
in Figures [l] and || (upper and lower thin solid lines), a significant change in 
the value of xf f leads to a significant variation in the resulting fractional error. 
Consequently, we have derived the value of that minimizes the maximum 
fractional difference in the whole mass range analyzed. The results, denoted by 
x^, are shown in column (6) of Table 1 for the case of S = 0.5, an intermediate 
value of the circularity parameter. The corresponding fractional error curve in 
Figures |l| and ^ has been plotted as a thick dashed line. The values of x^ are 
systematically smaller than the best fit values from S98 but always within the 
90% confidence intervals from Table 1. 



Table 1. Scaling of the NFW Density Law for Different Cosmogonies 



P(k) 


^0 


Ao 


0"8 


h 


m 
x sf 


X sf 


90% c.i. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


SCDM 


1.0 


0.0 


0.63 


0.5 


0.14 


0.17 


(0.12,0.25) 


ACDM 


0.25 


0.75 


1.3 


0.75 


0.18 


0.29 


(0.18,0.52) 


n = -1.0 


1.0 


0.0 


1.0 


0.5 


0.10 


0.17 


(0.10,0.30) 



* best fit value from S98 



Recent N-body simulations by Moore et al. (1997) suggest a profile of differ- 
ent shape from the one advocated by NFW. This could rise the question whether 
our results do depend on the form of the universal law assumed to describe dark- 
matter halo density profiles. To clarify this point we have repeated the previous 
analysis for different universal halo density laws in addition to the NFW pro- 
file. The latter corresponds to (a, (3, 7) = (1,1,3) in equation (0), while the 
rest of profiles considered are: the Hernquist (1990) law, (a, /3,t) = (1,1,4); 
the three-dimensional (3D) King law, (a,/3, 7) = (0,2,3); and the profile pro- 
posed by Moore et al. (1997) to fit their high resolution N-body simulations, 
(a, (3, 7) = (1.4,1.4,2.8). One can scale the values of x s and 5 C from the NFW 
profile to the rest of density laws by assuming that the maximum in the differ- 
ent circular velocity profiles is at the same physical distance to halo center. The 
Hernquist and NFW profiles are similar in the innermost regions of halos. The 
Hernquist profile has a different slope from the NFW profile in the outermost 
regions of halos although it was found to provide a satisfactory fit to halo density 
profiles (NFW). On the other hand, in the outermost regions the NFW and the 
3D King laws are similar. The 3D King law has a flat core which up to now 
has not been found in any cosmological N-body simulation (within the nominal 
resolution scale) , but this profile provides a good description of the distribution 
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Figure 3. The best predictions of our model for the SCDM cos- 
mogony in the case of the NFW density profile (solid line), the Hern- 
quist profile (dotted line), the Moore et al. profile (short-dashed line), 
and the 3D King law (long-dashed line). 
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of galaxies in clusters. For the present study, the King law can be considered as 
providing a limiting case for the slope of the density profile in the central region 
of halos. Since the exact value of this slope is still a subject of debate, we have 
included in our study the profile suggested by Moore et al. (1997). 

For each density profile we have derived the value of the dimensionless scale 
radius at formation that, according to our model, minimizes the maximum frac- 
tional difference between x e si and x* f in the mass range considered. As illustrated 
by Figure || corresponding to the SCDM cosmology, all density profiles lead to 
essentially the same result concerning the fractional difference. On the other 
hand, with respect to the value of inferred for the NFW profile, the values 
obtained for the other density laws analyzed differ in less than 20%. 

5. Conclusions 

The results of our study confirm the picture proposed in S98 and by Salvador- 
Sole et al. in these proceeding to describe the scaling evolution of universal 
halo density profiles. Halos form with a universal value of the dimensionless 
scale radius and scale density, x s f and 5 c f, i.e., halos have stable density profiles 
at formation. Subsequently, halos evolve by accretion and the corresponding 
evolution of x s and S c is well described by equations (|T|) and (|2|). Our study 
has covered a mass range of four decades, typically from the mass of galaxies to 
the mass of galaxy clusters. We have shown that the validity of the proposed 
picture is not affected by the specific form assumed for the universal density 
law of dark- matter halos. For density laws of the form (||) we always reach 
at essentially the same best model-predicted fractional difference between the 
values of x s f for halo progenitors and the new halos they form. We have derived 
for different cosmogonies the values of x s f which lead to the best consistency 
with the previous picture. The resulting values are found to be in agreement 
with the best fit ones obtained by S98 from the analysis of N-body simulations. 

According to the previous conclusions, in order to characterize the scaling 
evolution of halo density profiles, a first physical principle such as energy con- 
servation during halo formation has to be taken into account. But at the same 
time, one also requires a clustering model properly describing halo formation 
and evolution through an adequate distinction between major mergers, tracing 
the formation of new objects, and minor mergers, corresponding to accretion 
onto otherwise relaxed halos. On the other hand, the exact form of the uni- 
versal density profile of halos, which can not be determined by our analysis, 
would be set by more complicated processes than those included in the present 
study. Those processes would probably be linked to violent relaxation during 
halo formation in the framework of approximately self-similar cosmogonies. 
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